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Introduction 



It appears that, along with many of my friends and colleagues, I had been 
brainwashed by the great and tragic lives of Abel and Galois to believe that 
no general formulas are possible for roots of equations higher than quartic. 
This seemed to be confirmed by the brilliant and arduous solution of the 
general quintic by Hermite. Yet, below we find a formula giving a root to 
any algebraic equation of degree 2-5 and any reduced equation (see below) of 
^ ■ higher degree. This algorithm, which must have been familiar to Lagrange, 

^ . resulted when I was working on a paper on the asymptotics of hypergeometric 

functions where Gauss' multiplication formula for the gamma function is 
used to reduce certain infinite series, and by a happy accident my copy of 
Whittaker and Watson opened at p. 133. 



The Formula 

Without loss of generality it is sufficient to find at least one root to the 
reduced equation 

x^-x + t = {N = 2,3,4 ■■■)■ (1) 



^ ■ Letting x = ( ^^^^ ^\ we easily find that (1) becomes 

^' C = e'-^ + t0(C) (2) 

where 

Lagrange's theorem states that for any function / analytic in a neighborhood 
of a root of (2) 

/(C) = f{e'n + E -f^[/'(a)|0(a)r].=e-. (4) 

n=l ^- "''^ 



We now simply let /(C) = (^^^/(^^^), carry out the elementary differentiations 
(noting that D^x^ = T{p + l)xP~'^/T{p — k + 1)) and we come up with the 
root 

(N-2 further roots are found by replacing exp{27ri/{N — 1)) by the other N-1- 
st roots of unity, and the remaining root from the relation J2xj = 6^,2)- By 
the use of Gauss' multiplication theorem, the infinite series can be broken up 
into a (finite) sum of hypergeometric functions. 

_ -1 t ^V V^ ^^ )q„qN/(N-l) 1U=0 M JV ) 

' {N-iyUniN-i)^,^N-i^'' r(^ + i)n£o^r(2±^)^ 

giV/(iV -1) + 1 qN/iN-l) + N 
n+iFn[ , . . . , , 1; 

g + 2 q + N q tu .N-inrNi 

where u = exp{^^). In practice, n+i^n will always be reducible to at least 
nFn-1- Hence the root is a sum of at most A^ — 1 hypergeometric functions. 
The one technical point is that the convergence of these series requires that t 
be "sufficiently small", but this can be overcome by certain hypergeometric 
identities tantamount to analytic continuation. 



N=2 



Here we have 



Examples 



x^ -x + t = 



00 j.n 



However, by Gauss' formula 

r(2n + 1) = 4"(1/2)„(1)„ {{n)k = T{n + k)/T{n)) (7) 



so 

a;i = l-t 2^1(1/2,1; 2; 4t) (8) 

Since 

.F,(l/2.1;2;.)^?{;-^ j;j^; (8) 

we reproduce the quadratic formula. Note that the second root comes from 

a^i +a;2 = 1- 

N=3 

x^ -x + t = 

By separating the sum in (5) into sums over the even and odd values of 
n we obtain 

"^^ ~ ~ ^2^^r{n+ l)r(2n + 2) + 2" ^^^ r(n + 3/2)r(2n + 3) ' ^' 

By breaking up the gamma functions of multiple argument by using Gauss' 
multiplication theorem, the sums are easily identified as hypergeometric se- 
ries: 

t 3t^ 

xi = -1 - - 2^1(1/3, 2/3; 3/2; 27tV4) + — 3^2(5/6, 7/6, 1; 3/2, 2; 27t74). 
z o 

(10) 

However, from A. P. Prudnikov et.al. Integrals and Series, Vol.3 [Gordon and 

Breach, 1990]) we find 

2Fi(l/3, 2/3; 3/2; z) = ^ sin(i sin"! ^) 



iz 
18, .\ 



^F^ih/Q, 7/6, 1; 3/2, 2; z) = -[cos(- sin-^ v^) - 1], (11) 

Z o 

and we therefore have the three roots 

1 . .1 . _,.. r^,... .1 



xi = y= sin[- sin-i(tV27/2)] - cos[- sin-^(tv^/2)] 

V 3 3 3 

11 1 

X2 = y= sin[- sin-^(tV27/2)] + cos[- sin-^(tV27/2)] 

v3 3 3 



Xs = ^sm[^sm-\tVff/2)]. (12) 



Once again, for t > 2/v27 equation (10) must be analytically continued 
to obtain the correct form of (12). This amounts to writing sin"^ z = f — 
iLn{z + V-2^ — !)• 

Conclusion 

For N=2,3,4 Eq.(5) is definitely not preferable to the standard formulas, 
but for N=5, e.g. we get the root 

l/5,2/5,3/5,4/5;3125t^ 
" = '^^^f 1/2,3/4,5/4; ^5^^ ^'^^ 

in an elementary fashion with considerably less difficulty than by following 
the procedure in Davis' book [Introduction to Nonlinear Ordinary Differential 
Equations (Dover)]. It might also be pointed out that the above procedure 
carries over in a trivial way to the trinomial equation 

y^ - ay^-i + a = (14) 

where y = l/x, a = 1/t. Numerically, these formulas are not of much use 
since solutions can be obtained with the push of a button on many pocket 
calculators, but formulas such as (13) should have numerous entertaining 
uses, such as summing the odd hypergeometric series. 
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